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Abstract 




Interplanetary space provides simultaneously the best 
vacuum available to man and, because of the solar wind, a 
tenuous and unsteady high-speed outflow of predominantly 
hydrogen gas from the Sun, a remarkable variety of rarefied 
gasdynamics phenomena to observe. This paper provides a 
review of these phenomena, and of the way in which the pre- 
sent level of understanding has been achieved. 


Introduction 


Rarefied gasdynamic phenomena observed in interplane- 
tary space are not some insignificant stirrings of a few 
remaining particles in a near vacuum but a vital part p£ 
the transmission of effects of phenomena on the Sun to the 
Earth. It is not surprising that many aspects of inter- 
planetary gasdynamics were totally unexpected before the 
satellite era. On the other hand, certain features had 
been anticipated in theorizing about the way in which ob- 
servable solar and terrestrial phenomena might be related. 
From such studies in which observations were “explained*^ 
and "predictions" were made in terms of simple ^ ' 

are now well into a higher level of investigatiorf^h' which 
observations are more comprehensive and accura/^^and theor 
ries are more refined and quantitative. We it^^nd tete-Tfo ^ 
provide an account of how this level has beeru^^taihEc^ 
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some indications of the directions in which theory must be 
developed to keep pace with important observational develop 
ments. 


The Pre-Satellite Era 

For more than a century before the launching of the 
first spacecraft, it was known that a strong positive cor- 
relation exists between the sunspot number, auroral acti- 
vity, and small transient changes in the geomagnetic field. 
From this, it was evident that the Earth’s magnetic weather 
was influenced by changes on the Sun’s surface in contrast 
to ordinary weather which is influenced by the changing 
geometrical relationships of the Earth and Sun. 

In contrast to the slow and painstaking efforts 
required to establish statistical connections, the first 
indication of an explicit event on the Sun that might be 
connected with significant terrestrial consequences was 
sudden and unexpected. On September 1, 1859, the English 
solar astronomer, Carrington, saw an outburst of bright 
light within a large sunspot group; but within five minutes 
all trace of the event was gone, and the Sun appeared 
exactly as he had sketched it just before the event took 
place. At the same time, the Earth’s magnetic field became 
abruptly disturbed. About 18 hours later, a great geomag- 
netic storm unprecedented in intensity and duration com- 
menced abruptly. For several days, magnificent auroral 
displays were observed and telegraph communication was 
interrupted because of the current produced in the wires. 

In some cases this proved so powerful that the batteries 
were disconnected and the wires simply connected with the 
Earth. We would now regard the brief geomagnetic distur- 
bance on September 1 as due to the fleeting enhancement of 
ionospheric currents by ultraviolet light and X-rays; and 
the great storm that commenced the next day as the effect 
of an interplanetary shock wave followed some hours later 
by solar matter ejected fromthe Sun by Carrington’s flare. 

While Carrington cautiously proposed a connection be- 
tween the solar and terrestrial events, Kelvin and others 
dismissed the evidence as mere coincidence. Much later. 
Chapman and Bartels [1940] wrote in Geomagnet ism (a 1050 page 
two volume work published just one year after Chapman and 
Cowling’s Mathematical Theory of Nonun iform Gases, and , . 

lauded in reviews as the greatest on the Earth’s mag- 
netisra since Gilbert’s de Magnate published in 1600J) 
that Carrington’s observation of a flare in white light 
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remains unique, and that geophysicists would even then 
maintain a cautious attitude except for the clear relation 
found in 1935 and 1936 between solar eruptions and radio 
fade-outs which in turn tend to occur in coincidence with 
geomagnetic disturbances. 

From this uncertain beginning, the Norwegians, 
Birkeland and Stormer began in 1896 a prodigious, experi- 
mental, theoretical, and observational effort to explain 
the aurora and related geomagnetic phenomena in terms of 
beams of electrons emitted from the Sun and guided into the 
polar regions by the geomagnetic field. (See Stormer 
[1955) for a summary.) In 1911, Schuster showed, however, 
that electrostatic repulsion would disperse the beam to 
far too low a density to account for the terrestrial phe- 
nomena. Nearly a decade later, Linderaann proposed that 
geomagnetic storms might be produced by electrically neu- 
tral streams of charged particles; i.e., a plasma. 

Chapman soon thereafter was at work developing conse- 
quences of this idea, first alone and then in collabora- 
tion with V. C. A. Ferraro. 

By the beginning of the satellite era, there had 
emerged, (See Chapman [1963] for a rdsum^.) an incomplete, 
somewhat faulty, but still remarkably good, understanding 
of how ionized gas of solar origin might produce the ob- 
served geomagnetic variations and many associated effects. 
Interplanetary space was conceived to be primarily a va- 
cuum invaded at times by clouds of plasma ejected by solar 
flares, and possibly continuously by longlived beams of 
plasma rotating with the Sun. The rotating beams were used 
to account for the 27 -day recurrence pattern in geomagnetic 
and auroral activity by associating it with the rotation of 
the Sun. The plasma clouds were postulated to account for 
the flare-induced geomagnetic storms such as that associ- 
ated with Carrington’s observation. 

Chapman and Ferraro also established that the geomag- 
netic field would shield the Earth from the solar plasma 
and that a thin current sheet at the boundary of the plasma 
^^^Id terminate the geomagnetic field. They discussed the 
nature of the interaction with the front of either a flare- 
induced plasma cloud or a rotating beam, developed a pre- 
cise mathematical representation for both the steady and 
unsteady cases, and solved a number of problems to 
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illustrate particular features of the Interaction^ 

Near the Earth, the geomagnetic field is represented 
by a magnetic dipole at the center of the Earth having such 
a strength that | b| = = 0.312 gauss at the geomagnetic 

equator, and oriented so that the north geomagnetic pole is 
at 78.6^ North latitude and 70.1° West Longitude, near 
Thule, Greenland. Its properties are thus given by 


B 



(a^/r)^ 


yN A 

(0 sin 0 + r 2 cos 0) 


( 1 ) 


8 

in which a^ = 6,37 x 10 cm is the radius of the Earth, r 
is the geocentric distance, 0 is the polar angle measured 
with respect to the north geomagnetic pole and r and 0 
are unit victors in the r and 0 directions. 


No currents were presumed in the ^:sur rounding region, 
which is now called the magnetosphere^ hence £y£l § as 
well as div ^ vanish there. The condition chat the normal 
component Bn of the geomagnetic field must vanish at the 
boundary is supplemented by the relation B = p that 
equates the magnetic pressure at the boundary to the pres- 
sure p. of the incident plasma. Although the details of the 
interaction are very different from those in ordinary ra- 
refied gas flow incident on a solid obstacle, the net ef- 
fect in a plasma having an ion number density n and mass 
m flowing with a free-stream velocity v is that 
p - mnv^ cos'^V on an element of the boundary having its 
normal at an angle t to the free-stream direction. 

Although the model had been established much earlier, 
it was not until well into the satellite era that even 
approximate solutions were obtained for the shape of the 
geomagnetic field boundary [Spreiter and Briggs, 1962 and 
Briggs and Spreiter, 1963], Higher order solutions have 
been given most completely by Olsen [1969], but the results 
remain essentially the same as indicated by the earlier 
approximate solutions. In the idealized model of Chapman 
and Ferraro, the boundary is impermeable to the solar 
plasma, except, possibly at a pair of neutral points at 
which the magnetic field vanishes and from which extend the 
only field lines that connect the boundary and the Earth, 
These points are of considerable current significance 
because they define weak spots in the boundary through 
which interplanetary plasma can enter the magnetosphere 
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and precipitate into the upper atmosphere. 

In the years immediately preceding the launching of 
the first satellite^^ analysis of the behavior of comet 
tails by Biermann [1951, 1957] and theoretical studies by 
Parker [1958] began to suggest the continuous rather than 
intermittent existence of a high-speed flow of solar plasma 
throughout interplanetary space. The idea of such a solar 
wind did not gain acceptance, immediately, however, since 
the comet analysis involved a number of uncertain assump- 
tions, the theory was ambiguous, and no direct measurements 
could be made. 

The nature of the ambiguity in the theory can be un- 
derstood easily. The theory is just that of radially 
symmetric flow of a perfect gas in a centered inverse square 
gravitational field, under the assumption that the energy 
equation could be substituted for by assuming a simple 
polytropic law between the pressure and density. The 
governing equations have a family of solutions reminiscent 
of those for flow in a Laval nozzle. Parker chose the tran- 
sonic solution that takes the flow from low subsonic speeds 
near the Sun to supersonic speeds beyond a few solar radii. 
Chamberlain [1961] however, argued that purely subsonic 
solutions were the appropriate ones. A few years earlier, 
moreover, Bondi [1951] had used the same basic model, but 
with the direction of flow reversed, to explain accretion 
of mass by a star, and possibility the high 2 x 10^ 
temperature of the solar corona. 

Results of Direct Measurement of the Solar Wind 

It was Mariner 2 on its flight to Venus in the latter 
half of 1962 that resolved the controversies. The results 
[Neugebauer and Snyder, 1966, 1967] showed that the plasma 
flow was supersonic, always present, and that the mean 
values of the velocity, density, and temperature were 
approximately as indicated by the ®'first order” solar wind 
theories then available. With the existence of the solar 
wind established, there followed many refinements of the 
theory in which effects of heat conduction, viscosity, 
rotation, magnetic fields, and, most recently, fluctuations 
of the flow were included. 

Later and more complete measurements with instruments 
that scanned in several directions have shown that the 
velocity distribution of the random motions of the particles 
is not isotropic and that the ^’temperature” is therefore 
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different in various directions^ By defining a quantity 
analogous to ordinary temperature in terms Of the random 
motion of the particles in a given direction and by refer- 
ring to it as the temperature in that direction, Hundhausen 
et al. [1967] have shown that the direction of the tempera- 
ture anisotrophy is characteristically aligned with the 
interplanetary magnetic field, the maximijm temperature is in 
the direction pointing away from the sun along the field 
lines, and the minimum temperature is in the direction trans 
verse to the magnetic field* This indicates that heat is 
being conducted away from fae sun along the magnetic lines of 
force* It is evident that such anisotropies have an 
important effect on the dynamics of the solar wind, but 
relatively little consideration has yet been given to them 
in refinements of the theory* 

Another feature that needs additional consideration is 
the rotating sector structure discovered by Wilcox and Ness 
[1965] in the IMP 1 data* Although the sector is defined in 
terms of the predominant direction of the interplanetary 
field as away from or toward the sun, they showed that the 
velocity, density, magnetic field intensity, and geomagnetic 
activity index all display an organized reproducible varia- 
tion within each sector* Wilcox [1968] and Schatten [1971] 
have examined the evolution of the sector structure through- 
out the years since its discovery* From data observed over 
less than whole 11-year solar cycle, and inferences from 
ground-based geomagnetic data, they conclude that the four 
sectors found in the IMP-1 data are representative of the 
solar minimum, whereas two sectors are more representative 
of the period near solar maximum* 

In addition to the large-scale irregularities associ- 
ated with the sector structure, the solar wind exhibits a 
wide variety of variations, both gradual and abrupt. Many 
are simply convected with the flow* Others propagate 
through the solajx wind as shock waves* Most of these may 
be identified with large solar flares that occurred a few 
days previously. However, most solar flares do not produce 
shock waves at the orbit of Earth* Hundhausen [1972] used 
data from 19 shock waves, to seek the cause for such a 
selection* He concluded that the observed shock waves are 
produced by the relatively few flares that have a suffi- 
ciently high energy-mass ratio for escape against solar 
gravity* 
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Magnetohydrodynamic Representation of Solar-Wind Flow Fields 


Most analyses of structures in the solar wind, and of 
the flow fields of the Earth, Moon, and planets, have been 
based on the continuum equations of gasdynamics or magneto- 
hydrodynamics of a perfect dissipationless gas. That they 
can describe the behavior of a gas so rarefied as the solar 
wind is one of the surprises of space research, but their 
appropriateness has been confirmed in a wide variety of 
applications. To be precise, the equations are as follows 

dp/9t + ^ . pv = 0 (2) 

pDv/Dt = - VP - B X curl B/4 tt + pg (3) 


SB/^t = curl (v X B) , 


Ds 

Dt 


0 , s-s 

o 


c In 

V 


div B = 0 


p/Pq 

(p/Po>^ 


(4) 

(5) 


in which p, v, B, S, and g represent the density, velocity, 
magnetic field, ^entropy, and gravitational acceleration. 

If effects of gravity are disregarded, solutions depend on 
two nondimens Lonal ratios, Mach number M = v/a and Alfven 
Mach number = v/A, where a = (vp/p)^/2 is the speed of 
sound and A = (B^/4np)^'2 is the Alfven speed. Both M and 
are normally much greater than unity, values of about 10 
being representative for the solar wind flow approaching 
the Earth bow wave. 


These equations must be supplemented by additional re- 
lations between conditions on opposite sides of possible 
discontinuities in the flow. At each element of such a sur^ 
face, conservation of mass, momentum, energy, and magnetic 
field provide 

Ipv^J = 0 , [pv^v + (p -h B /8fr)n - B^^/4rr] ^ 0 (6) 
[pv^{h + V /2) + v B /4rr - (B v « B)/4n] = 0 (7) 

^ A. 

'Vc - ?t ”nl ■ 0 . ■ 0 (8) 
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The square brackets indicate the difference between the en- 
closed quantities on the two sides of the discontinuity; 

V = Y - n - ^ is the normal fluid velocity component re- 
lative So the normal velocity \ of the discontinuity sur- 
face; n and ^ are unit vectors normal and tangential to 
the discontinuity surface; a.nd subscripts n and t indi- 
cate components of v and g in these directions. Five 
classes of magnet ohydrodynaraic discontinuities, tangential, 
contact, rotational, and fast and slow shock waves, are de- 
scribed by these equations. All except the contact discon- 
tinuity have been identified in the solar wind, and several 
are important. in the flow field of the Earth and other ma- 
jor objects in the solar system [Hundhausen, 1970, and 
Spreiter and Alkane, 1970]* 

An important refinement of this theory is currently 
being developed to explain a number of .features of the 
solar wind associated with the observea anisotropy of the 
temperature and pressure. In this approximation (see 
Burlaga [1971] for a review), the density p and velocity 

V are defined by o = E m.n, and v = i L m.n.u, in which 

^ ^ LL ,«pLl~L 

m. is the particle mass; and n^ , the number density of 
particles of species i, and the average speed of spe- 
cies L, are given in terms of the velocity distribution f 
in the usual way by n(r,t) = /f(r,v,t) dv, u(r,t) = i 

fy, M A# n 

/ vfdv . The governing differential equations remain as 


written in equations (2) through (5), except that the right- 
hand side of the latter is replaced by 

V.[I(2 P, + B^/8n) - B B |/4n] + pg 
~ k ~ ~ 

which reduces to 


k 1 


B^/Bn) 


- (B.7) B 5/4tt + pg 


^^a measure of the anisotropy, 

4i 


when 

ter quantity is defined as g 


= 1 - 


Pk„ ' "A I ”k„”k„ ' 


is constant, 

3 ^ ’Pk 

B k 1 


The Lat- 
in which 


and Pj^ = Wj^ f dv 

-L th 

represent the pressure of the k species parallel and 
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perpendicular to B , and w, = v, - u, . These equations 

A# A» 

imply collective interaction of the particles# For motions 
normal a-o B, this is caused by the magnetic field; colli- 
sions are not necessary. For motions along B, details 
are not yet fully understood, but it is believed that such 
interactions are associated with fluctuations or wave-par- 
ticle interactions# It is possible, moreover, that on 
Some scales the interactions along the field lines are so 
weak that one should use a kinetic particle equation for 
motions along B. 

Equations (2) through (4), or their anisotropic coun- 
terparts, must be supplemented to provide a closed set of 
equations# For the isotropic case, only one additional 
equation is needed, and it is frequently sufficient to use 
the adiabatic approximation indicated by equation (5)# 

For the anisotropic case, there is one more dependent 
variable because of the two pressures pn and p , and it 
is necessary to supply two additional reiationsf 

One approximation that has sometimes been applied to 
the solar wind is that of Chew, et alo [1956] in which 
^(P /pS)Dt = 0 and D(piiB^/p^) /Dt - 0 # The first is sug- 
gested by an adiabatic* 'invariant of charged particle mo- 
tion in a magnetic field; the second stems from the as- 
sumption of either no or constant heat flux along JB# 

These equations were derived on the assumption that g =? 
nkT/(B /8tt) « a condition that seldom applies in the 
solar wind. 

For applications to the solar wind, Burlaga [1971] 
proposes chat p^ = pYj^ and p^ ' ^ P^ll which ap- 

II II 

propriate values for yn and y are still to be deter- 
mined. For relatively ‘'small-scale features of the solar 
wind for which the magnetic field is reasonable orderly, 
he proposes that y =2 because there are 2 degrees of 
freedom normal to B and yjj = 3 because there is only one 

degree of freedom along B- If, however, compressions along 
B are isothermal, since the conductivity along B is high,' 
fie suggests that perhaps y.i might be 1 rather than 3. It 
is quite possible, furthermore, that the appropriate value 
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for Y depends on the scale of the phenomena being consid- 
ered, since the magnetic field is probably more disordered 
when viewed on a large scale than on a small scale. 

Just as for the isotropic fluid, additional equations 
must be supplied to relate conditions on opposite sides 
of discontinuity surfaces. We will not present further 
details, but refer the reader to Burlaga [1971], There 
may be found, in addition, an account of the properties of 
both waves and discontinuity surfaces in the anisotropic 
collisionless medium described above, and a discussion of 
a variety of observations in space that give a sense of 
reality to the theory. Although the theory is clearly 
applicable to the solar wind, many consequences remain to 
be worked out. Nonlinear theories of waves, instabilities, 
and interaction of discontinuities aqd waves are needed, as 
are further developments of the effects of fluctuations, or 
possibly turbulence, in the flow. The effects of anisotro- 
pies on large-scale features of the solar wind, such as the 
sector structure, and on the flow about the Earth and other 
major objects in the Solar System are almost totally un- 
known at present, and demand investigation if one is to 
have a consistent theory of significant solar wind pheno- 
mena. 

Solar Wind Flow Past the Earth - A Magnetic Planet 

To calculate solar wind flow past the Earth, a know- 
ledge of the density and velocity of the solar wind is 
needed. Even before this was acquired by Mariner 2 in 1962, 
early satellite data had confirmed the termination of the 
geomagnetic field at about the expected location. They 
also revealed an unanticipated transition region having 
fluctuating properties before steadier conditions typical 
of interplanetary space were reached. Various hypotheses 
were put forth at the time; but Axford [1962] and Kellogg 
[1962] correctly explained that the Earth’s magnetosphere 
had a bow shock wave. Shortly thereafter, Spreiter and 
Jones [1963] used gasdynamic theory to calculate its loca- 
tion, They also noted that the fluctuating character of 
the post -shock gas Is not unlike that calculated for 
©Dllisionless shock waves by Auer et alo [1961, 1962], 
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The first Earth satellite to repeatedly carry plasma 
probes and magnetometers across the geomagnetic boundary 
and the b.ow wave was Explorer 18 or IMP-1 launched on 
ovember 27, 1963, As in the earlier measurements, the 
raasLtion region is recognizable in the magnetometer re*> 
cords by the presence of fluctuations more substantial 
^ incident solar wind or the magnet os phereL 

e differences are even more obvious in the data from the 
plasma probes. aU [1966] shoved there was a to- 

ta absence of plasma flux when the spacecraft was within 
the magnetosphere and shielded from the flowing solar plas- 
ma. The presence of plasma flux in only a narrow range of 
energy levels satellite was beyond the bow wave indicates 
that the random thermal velocities of the particles in the 
incident solar wind are small compared with the directed 
bulk velocity of the flow. The broad energy spectrum ob- 
served in the intervening part of the orbit is indicative 
of the hot shocked plasma in the transition region. 

Although the magnetopause and bow wave were not al- 
ways distinctly revealed it was usually possible to dis- 
tinguish the various regions* They were, moreover, close 
to the theoretical positions of the magnetopause and bow 
wave calculated by Spreiter and Briggs [1962] and Spreiter 
and Jones [1963], Although the theoretical model was an 
inconslstant mixture of particle and fluid concepts, the 
results served to fix the ideas about the fluid-Uke nature 
of the flow field, and to encourage a re-examination of the 
entire problem from a consistant magnetohydrodynamic point 
of view. 

The problem of solving the magnetohydrodynamic equa- 
tions (2) through (8) for steady uniform flow at infinity 
past a magnetic dipole described by equation (1) remains 
intractable j but Spreiter et al. [1966] showed that approx- 
imate Solutions of good accuracy can be attained upon in- 
troduction of several important simplifications beyond the 
obvious elimination of terms containing S/5t or X , 
Within the magnetosphere, B^/grr greatly exceeds p every- 
where above a few hundred kilometers; and the dominant ef- 
fect of the Earth is provided by the terms that remain wher 
p and p are equated to zero, namely div B = 0 and curl B 
0< Near the Earth, B may be represented adequately by 
equation (1), The magnetopause must be represented by a 
tangential discontinuity, since that is the only solution 
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of the conservation equations for which = 0* Ar- 
bitrary differences in p ^ , and are allowed across 

such a surface, but p + must be the same on both 

sides. The bow wave can be represented only by a fast- 
shockwave solution of the conservation equations because 

the solar wind approaches the Earth with a mass flux p v 

■^00 00 

that greatly exceeds that which can pass through any of the 
ocher types of discontinuitieso With the neglect of the 
small effects of terms containing B and g in equations (3), 
(6), and (7) , the equations for the fluiS motion reduce to 
those of gasdynamics; and the distortion of the interplane- 
tary magnetic field can be determined in a subsequent step 
by solving the remaining equation with v known from the 
gasdynamics results. 


The free-boundary problem for the shape of the magne- 
topause can be made tractable by introducing the Newtonian 
approximation for the pressure on the magnetopause, i.e. by 
P ^ = B /Srr in which, K is a constant equal 

to 0.88 for high Mach number flow of a monatomic gas, al- 
though usually taken as unity in most applications of type 
described here. Since the resulting description coincides 
with that put forward more than 40 years ago by Chapman and 
Ferraro., it is paradoxical, but true, that the theory of 
the geomagnetic boundary shape is much older than the idea 
of the solar wind. It follows that the locations of the 
magnetopause and bow wave calculated previously remain ap- 
propriate, but the logical inconsistencies have been re- 
moved from the underlying theory. 

The most detailed comparison of results calculated in 
this way and direct observations in space are those made by 
Spreiter and Alksne [1968] utilizing data from Pioneer 6 
following launch on December 16, 1965, during a period of 
exceptionally low geomagnetic activity. The calculated and 
measured locations of the magnetopause and the bow wave 
were shown to be virtually coincident, and all the flow pa- 
rameters to be in remarkable agreement with the theory. It 
is difficult to assess the implications of the remaining 
discrepancies, however, because of uncertainties and incom- 
pletness of the data. The experimental values for p, for 
example, were judged by the experimenters [Wolfe and 
McKibbin, 1968] to contain a 50 percent uncertainly, and 
these given for the velocity are not actually for the bulk ^ 
velocity, as considered in the theory, but for the speed of 
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the Lons that produce the greatest current per unit energy 
increment in the plasma prnbeo It is possible, moreover, 
that some of the observational variationa not duplicated 
by the theory are the result of unkown changes in the in- 
terplanetary conditions that occurred when Pioneer 6 was 
behind the bow wave, and could be reproduced by the theory 
if the actual interplanetary conditions were better known. 

At the present time, we may summarize the present 
state of the magnetohydrodynamic theory of the interaction 
of the Solar wind and the Earth as being most thoroughly 
worked out for the more upstream portion of the flow field 
and for solar wind conditions that are steady and charac- 
terized by higti (5 or greater) values for both MAc# to ^ • 
In addition to the obvious need for solutions for the dyna- 
mical response of the flow to abrupt changes in conditions 
in the incident stream, since these often produce the most 
dramatic geophysical consequences, there are many ways in 
which the theory is notably incomplete. First of all, low 
values for , occasionally less than unity, are some- 
times observed in the solar wind. There only now is begins 
ning a discussion of such possible effects on the flow about 
the Earth and the attendant geophyical consequences 
[Formisano et al, , 1971; Rizzi, 1971; and Fairfield, 1971] - 
There is a great "need for improved theoretical understand- 
iiag of the properties of the magnetic neutral sheet and 
the enveloping plasma sheet that extends downstream along 
approximately the centerplane of the magnetosphere tail. 

As noted above, the aeutral points on the magnetosphere 
boundary have been long suspected from suggestions of 
Chapman and Ferraro to be of importance as exceptional re- 
gions through which charged particles from the solar wind 
could gain access to the inner magnetosphere. Spreiter and 
Summers [1967] have considered implications of the fluid 
theory for this region, concluded that a cusped shaped re- 
gion of hot "stagnant** plasma extends toward the Earth from 
the vicinity of the theoretical neutral points of the ide- 
alized theory, and given expressions for the leakage rate 
of charged particles from the cusp based on Grad's [1963] 
analysis of the particle leakage from the ends of a mirror 
machine. At present time, there is a considerable effort 
in progress to measure particle fluxes from these regions, 
and to ascertain the resultant geomagnetic, auroral, and 
ionospheric effects. The stability of the magnetosphere 
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boundary, and the bow wave,has been the subject of consider- 
able speculation, but a comprehensive theory is lacking, 
as it also is for effects of fluctuations and anisotropies 
of the flow„ Although much has been accomplished, the im« 
portance of acquiring a better understanding of the varied 
and subtle ways in which solar and terrestrial phenomena 
are related, demands that more be doneo 

Solar Wind Flow Past Mars and Venus - Nonmagnetic Planets 

If we now turn attention to Mars and Venus, we find 
that data acquired in recent years disclose a somewhat 
different type of interaction from that just described for 
the Earth, Neither of these planets has a significant mag- 
netic field to withhold the solar wind, but they both have 
an upper ionosphere that is sufficiently conducting elec- 
trically to prevent the solar wind from flowing directly 
into the lower absorbing levels of? the atmosphere. The 
solar wind is thus deflected around the ionosphere, and a 
bow wave is formed upstream of the planet, similar in many 
ways to that associated with the Earth, Aside from evident 
differences in the underlying physical processes at the 
ionopause, the surface that bounds the ionosphere and the 
Solar wind, the principal difference between the flow 
fields around these planets and the Earth is the size of 
the cavity. Compared to the Earth, for which the nose of 
the magnetopause is at an altitude of about 60,000 km for 
representative solar wind conditions the nose of the iono- 
pause is at an altitude of about 500 km for Venus and only 
about 155 to 175 km for Mars, 

It has been shown [Spreiter et al,, 1970] that a 
theory for solar wind flow past Venus or Mars can be con- 
structed analogous to that for Earth on the basis of the 
facts that the planetary magnetic field is weak or nonexis- 
tent, and that the ionospheric pressure p is sufficient to - 
stop the solar wind. If effects of gas motions are ne- 
glected, conditions in the ionospheres of these planets may 
be Idealized to the hydrostatic support relation grad p = 
p| , or dp/dr -pg that remains when v and B are'^zero in 
equation (3).,. If the scale height H, ="'kT/mg"^of the upper . 
atmosphere is constant, the tangential discontinuity 
becomes p = Kp^ cos^^ = exp [-(r - T|^)/H] 
in which is the pressure at a reference radius « 
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Once the ordinary differential equation for the shape of 
the ionopause is solved, the properties of the surrounding 
flow may be calculated in the same way as for the Earth „ 

Spreiter and Rizzi [1972] have compared the locations 
of the ionopause and bow wave calculated in this way for 
conditions representative of Mars with the magnetic field 
measurements of Mariner 4 [Smith, 1969] as it flew by that 
plane to They showed that Mariner 4 crossed the theoretical 
location of the bow wave at two points, and that the mag- 
netometer record displayed prominent discontinuities at 
virtually the precise times that the shock crossings were 
indicated. While there is always the possiblity that such 
sudden changes in the magnetic field may be the result of 
the passage of discontinuities in the solar wind, the coin- 
cidence with the theoretical results supports the proposi- 
tion that these data do indeed provide the signature of a 
Martian bow wave. A half-hour gap in the data during the 
two hour Mariner 4 was behind the bow wave is indicative of 
the many difficulties of space research. During this inter*- 
val, the spacecraft was directly in line behind the planet, 
as viewed from Earth, and no data could be received from it^ 
Moreover Mariner 4 carried a plasma probe, but a malfunc- 
tion prevented It from providing an independent indication 
of the shock crossings. Mariners 6 and 7 subsequently ap- 
proached within 2000km of the Martian surface, but neither 
spacecraft carried a magnetometer, plasma probe, or ener- 
getic particle detectors. Had the geometry of the Mariner 
4 encounter been slightly different, we might still have no 
direct observational knowledge of the nature of the Martian 
interaction with the solar wind and of the insignificance 
of its magnetic field in comparison with that of the Earth 

Spreiter et al., [1970] have also made a similar 
comparison for Venus using the data from Mariner 5. 

[ 1967 ] 0 They showed that Mariner 5 crossed the bow wave 
at very nearly its calculated location but that there is 
no positive indication the spacecraft penetrated the iono- 
pause to the extent indicated by the theory. Since the 
observations near the theoretical location of the iono- 
pause Indicate values for n/n<x> and v/v^ that are about 
one-half the theoretical values for the flow exterior to 
the ionosphere, Spreiter et al., [1970] suggested that per- 
haps Mariner 5 entered a thick boundary layer separating 


511 



JOHN R, SPREITER AND ARTHUR W. RIZZI 

““ '»- 10"o- 

Venus should have a thick boundary layer alone 
iahi whereas the Earth's magnetopause -is aston- 
diffaSnf plausible in view of the fundamentally 

rent nature of the two boundaries « The magnetopause 

rrelac“"^*^^^ ® boundary between the flowing plasma and 
tween the ionopause is ‘a boundary be- 

cation of plasma in relative motion. At the lo- 

tSn?La^J ^o^opause indicated by the present dissipa- 
hi plasma velocity might be expected '^to 

iLeeH ^ indicated by the theory, as is 

indeed evident in the data of Mariner 5. A more detailed 

b,“ ^ 

nient can be made, however* 

aminatL^nr carried out a more complete ex- 

amination of one aspect of the interaction by developing an 

’^he interaction unS 

nL dSectJo^ Sn^Lf ® aligned with the 

?■ incident solar wind, as it approxi- 

mately was at the time of the Mariner 5 encounter witT 
Venus [see Riazi (1971) for a complete account], so 
doing, the ionopause was still represented by a tangential 
discontinuity, but the calculated values for p + b2/|tt of 

Ld tie'florf^nt^ Newtonian approximation, 

computed using the equations of 
raagnetodydrodynamics rather than gasdynamics. 

These calculations were made possible by an extension 

forfat”^ embedded shock waves of a trans- 

formation scheme described for a succession of increasing- 
ly complex magnetodydrodynamic cases by Cowley (19601 ^ 

rSi' 1 7 1 L 

Rizzi [1971]. According to this procedure, a general pro- 

TToLli t f ~ ' ^hlch \(t) is 

a constant along each streamline, is supplemented bv a set 

p* . p * bX. b. Jb-. 

^ which p/4tt — ^ p/p . Substitution 

of these relations into the equations for steady magneto- 
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hydrodynamic flow yields equations for the transformed var- 
iables that are identical to equations (2), (3), (5), (6), 
and (7) with \ , and B equated to zero. They need 

to be supplemented by an ‘’equation of state*' relating , 
p* , and s* to complete the system. Although the equation 
of state relating p , p , and s is the same as in ordinary 
gasdynamics for a perfect gas, the corresponding relation 
between p* , , and s* is substantially different. It is 

that 


P* (p,s*) - (1 - p 


Y"1 p ^ ^co 


+ h 


CO A 


( — ) exp 
Pco 


in which p = p*/(l+X^p*/4n) = P*/( 

1 2 1 *2 ® 
= *7V + h =“:rv +h* = h* = cT 

^ CO Z OD OD St P CO 


s-s 



V 


and h 

St 


In this way, a correspondence is established between 
equations for magnetohydrodynamics for aligned flow and 
those of gasdynamics for a hypothetical gas obeying an un- 
usual equation of state. The advantage is that existing 
methods for solving the gasdynamic equations can be applied 
with only minor modifications to obtain magnetohydrodynamic 
solutions. 

Comparison of the results calculated in this way with 
the observations of Mariner 5 near Venus show that the 
shock location is predicted almost perfectly and 
that the penetration of the trajectory is not quite so deep 
into the inosphere* The conclusion is that a distinct im- 
provement in the quality of the theoretical representation 
has resulted from an obvious, but complicated, improvement 
in the theory. 


Solar Wind Flow Past the Moon 

Solar wind flow past the Moon is notably different 
from that for either the Earth or the nonmagnetic planets. 
Mars and Venus. As with the other objects, it was only 
with the acquisition of data measured in situ that an end 
could be put to a period of speculation about the nature 
of the interactiono Although more recent measurements are 
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revealing a rather rich substructure, plasma probe and mag- 
netometer data from Explorer 35 placed into orbit about the 
Moon in July 1967, provided immediately the knowledge need- 
ed to specify appropriate bouhdary conditions to complete 
an idealized magnet ohydrodynamic representation of the in- 
teraction. These conditions are that free«stream condi- 
tions exist directly up to the lunar surface, that the 
magnetic field is continuous from the solar wind into the 
Mo on^ that no significant electrical currents flow in the 
Moon in the steady state, and that there exists a void in 
the solar wind downstream from the Moon in which neither 
particles nor electrical currents are to be found. Two al- 
ternatives exist for describing conditions at the surface 
of this void, depending on whether or not an electric cur- 
rent sheet forms. If one does form, the boundary of the 
void must be represented by a tangential discontinuity. If 
no current sheet forms, the magnetic iiie Id in the Moon and 
the trailing void must join continuously with that in the 
surrounding flow. Both theory and data from Explorer 35 
have been shown [Spreiter et al, 1970] to indicate that 
either possibility may occur, depending on the orientation 
□f the interplanetary magnetic field. 

Both a discrete particle representation (Whang, 1968a, 
b] and a continuum magnetohydrodynamic fluid [Spreiter et 
al, 1970] have been developed .in detail to account for many 
aspects of the interaction. The latter reference, and also 
Ness [1970], provide an extensive discussion of the obser- 
vations of Explorer 35, both alone and in conjunction with 
the data from Explorer 33 which was simultaneously monitor- 
ing the solar wind away from the influence of the Earth and 
Moon, and their relation to the theoretical predictions. 
Small perturbations of the interplanetary magnetic field 
Intensity of the order of 20 to 30 percent are observed to 
be correlated with the location of the solar wind plasma 
umbra and penumbra* On some tranverses, characteristic al- 
ternations in sign (+-+-+) of the magnetic * perturba- 
tions are observed as the satellite traverses the lunar 
wake. On other traverses, no' such alternations are evident, 
and only slight variations are apparent. , 

Some of the variations are understandable in terms of 
the direction of the interplanetary magnetic field and the 
orientation of the orbit with respect to the principal 
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plane of symmetry defined by the incident solar wind velo- 
city and magnetic field directions [Spreiter et al. , 1970]. 
Others may be correlated with the diamagnetic properties of 
the solar wind, as measured by the ratio of the components ' 
of the plasma and magnetic pressures iperpendicular to the 
magnetic field direction [Ness^ 1970]. 

It is becoming evident from the data from the Apollo 
magnetometers on the surface of the Moon that the story is 
still far from complete, however. As reported by Driscoll 
[1972], the magnetic field at the lunar surface was mea- 
sured to be 38 gammas at the Apollo 12 site, 43 and 103 
gammas at two locations near the Apollo 14 site, only 6 
gammas at the Apollo 15 site, and 120, 125, 180, 230, and 
313 gammas at several locations near the Apollo 16 site. 
Since such intensities greatly exceed the critical value 
B = (STTKp^v 2) 1/2 required to stop the solar wind at a 
stagnation point, it is apparent that modifications may 
have to be made to all the theories that assume that the 
solar wind flows without deviation into the lunar surface. 
Since even weaker magnetic fields would be sufficient to 
deflect the solar wind away from the flanks of the Moon, as 
seen from the Sun, it is quite possible that many of the 
small increases in magnetic field intensity seen just out- 
side the lunar wake may be the result of localized outward 
deflections of the flow and their accompanying bow waves. 

It is apparent that future studies of the solar-wind lunar 
interaction will have to take these local magnetic irregu- 
larities into account- 


while it is true that all of the perturbations of the 
solar wind resulting from interaction with the Moon are 
small, their implications go substantially beyond being 
just one more source of disturbance in an already fluc- 
tuating solar wind. From these studies can come many in- 
ferences about the lunar interior, a matter of interest 
both for its own sake and also for the perspective it of- 
fers on our knowledge of the Earth’s interior. . As summa- 
rized recently by Ness [1970], for example, it has been 
concluded from these data that neither the Gold-Tozer-Wil- 
son mechanism of accretion of interplanetary field lines 
or the Sdnett-Colburn-Hollweg mechanism of unipolar induc- 
tion is of significance for the Moon. From the fact that 
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there is little distortion in the transmission of micro- 
structural discontinuities In the solar wind flow past the 
Moon, it has been concluded -that the lunar interior must 
have a low electrical. conductivity ■ (less than lO**^ mho/m) 
and ^erefore a relatively cool interior of the order of 
1000^, Reconciliation of such conclusions with the pre*- 
eence of surface magnetic field intensities of the order 
measured by the Apollo magnetometers, and with the appar- 
ently universally, high remnant magnetization of lunar 
rocks, is not easy; and in fact poses a major problem for 
anyone attempting to construct an acceptable model for the 
structure, composition, and evolution of the lunar inte- 
rior. It is an intriguing thought for this Symposium that 
studies of rarefied gasdynamics might contribute to under- 
standing the state of the lunar interior, not only at the 
present time but perhaps a few billion years ago. 

Solar Wind Interaction with Other Major Objects in the 

Solar System 

At the present time, direct measurements have been 
made of the interaction of the solar wind with the Earth, 
the Moon, Venus, and Mars in descending order of 
completeness. Although the surprises encountered with 
each of these objects should be sufficient to give pause 
to anyone about to speculate on the nature of the interac- 
tion with other major objects in the Solar System, the 
knowledge gained from their study gives a sense of perspec- 
tive about the possibilities that did not exist as recently 
as only a decade ago. Briefly, the outlook appears as fol- 
lows. 

Mercury has no significant atmosphere, and its slow 
rotation rate of 59 days suggests that it may have no mag- 
netic field. Its interaction with the solar wind is anti- 
cipated, therefore, to resemble that for the Moon, Jupiter 
emits radio signals that are interpreted as indicating that 
planet has an enormous magnetic field with a dipole moment 
that exceeds, that of the Earth by a factor of about 3 x 10 • 

It Is anticipated, therefore, to have a magnetosphere like 
the Earth, but scaled up in size. by a factor of approxi- 
mate ly 50, considering both its magnetic field and the an-^ ■ 
ticipated variation of the solar wind properties between 
the orbits of Earth and Jupiter. If, however, the solar 
wind should become subsonic inside the orbit of Jupiter, 
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there would, of course, be no bow wave, and a new type of 
Interaction would be antlcipated« Virtually nothing can 
be said about the interaction of the solar wind with the 
remainder of the planets because of almost total lack of 
knowledge of properties of both the distant solar wind and 
the planetary magnetic fields and ionospheres. An addi- 
tional possibility that, might occur in the outer part of 
the Solar System, particularly if the solar wind should be- 
come subsonic, is that the flow may be deviated primarily 
by the gravitational field. In such a case, the object 
would act rather like a sink in the solar wind, drawing in 
the plasma as it approaches. 

Noteworthy also is that five planetary satellites, in- 
cluding the four large Gal il lean moons, two of which are 
approximately the size of the Moon, are deeply embedded in 
the presumed magnetosphere of Jupiter, whether the solar 
wind is supersonic of subsonic. Similarly, our Moon spends 
about 7 days out of every 29.5-day synodic period in the 
region disturbed by the Earth. It is clear that conditions 
surrounding the Moon during these times differ greatly 
from those in the undisturbed solar-wind flow, but we are 
aware of no theoretical analysis of the details of the in- 
teraction, . 

. Yet another class of Interactions with the solar wind 
is that with comets having Type I, or ionized, tails. As 
noted at the beginning of this paper, it was, in fact the 
study of motions of such cometary tails that led to the 
first suggestion of the permanent existence of the solar 
wind by Biermann. While direct probing of a comet has not 
yet been attempted with a spacecraft, the great visibility 
of comets make them suitable objects for optical observa- 
tion from the Earth. Although the analysis is still a ra- 
ther preliminary stage of development, a magnetohydrodyna- 
mic model is beginning to emerge that is very similar to 
that described above for the Earth and neighboring planets, 
except that the cometary nucleus of 1 to 10 km radius is 
considered to be the source of a flow of neutral molecules 
as it comes into the warmer inner parts of the Solar Sys- 
tem. These do not interact with the solar wind until they 
are ionized at some distance from the nucleus primarily ' 
by photoionization and charge exchange, but probably by 
other mechanisms as well. The effect is similar in many 
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respects to that of distributed mass injection into the 
flow surrounding the comet, and leads to an interaction re- 
gion that is enormously larger than that which would be in- 
dicated by simple estimates based on solutions for flows 
without mass, injection- - 

Concluding Remarks : 

In conclusion, perhaps the outstanding impression to 
be gained from this review is a realization of the tremen- 
dous gain in understanding of the properties of the solar 
wind and its interaction with objects in the Solar System 
that has been achieved in recent years through an intensive 
observational and theoretical program- The best progress 
has been made when at least a minimum of observational evi- 
dence is in hand to guide the theory along realistic paths, 
and when there is enough theory available to provide rea- 
sonably correct interpretations of data that are character- 
istically incomplete in coverage and often imprecise in de- 
tail. Much has been learned about how effects of events on 
the Sun are conveyed to Earth by the interplanetary medium; 
and much has been learned about the dynamical behavior of 
rarefied ionized, gases. 6tudy ot the discrepancies between 
theory and observation is providing the inspiration to ex-; 
tend the theoretical analyses to include many features, 
such as those ' that :,involve effects of anisotropies, insta-, 
bilities, waves, discontinuities, and other fluctuations 
possible in an ionized rarefied conducting medium, that are 
beyond the reach of the current theories based on magneto-^ 
hydrodynamics or gas dynamics- Multicomponent collision- 
less plasma theory holds promise of providing understanding 
of many of these, but the required theory is much more com- 
plicated and far less developed for relevant flow configu- 
rations than the fluid theories reviewed here. Finally, 
perhaps the most intriguing fundamental question relates to 
the underlying reasons for the success of the continuum 
theories in dealing with the most rarefied gas presently 
accessible to direct measurements.. Its proper resolution 
must be ranked among the more: important goals to be achieved 
if the environment provided by interplanetary space is to . 
fulfill its often expressed potential as a new laboratory 
for Che study o,f rarefied plasmas- We thus conclude with 
the observation , Chat the study of interplanetary dynamical 
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processes provides both a new field of application for ra 
efied gasdynamicB, and the opportunity to obtain deeper 
insight into the foundations of the subject itself. 
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